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ABSTRACT 

We derive an analytic formula of genus for dark matter halos assuming that the 
primordial mass density field obeys the random-Gaussian statistics. In particular, we 
for the first time take account of the nonlinear nature of the halo biasing, in addition to 
the nonlinear gravitational evolution of the underlying mass distribution. We examine 
in detail the extent to which the predicted genus for halos depends on the redshift, 
smoothing scale of the density field, and the mass of those halos in representative cold 
dark matter models. In addition to the full model predictions, we derive an explicit 
perturbation formula for the halo genus which can be applied to a wider class of biasing 
model for dark matter halos. With the empirical relation between the halo mass and 
the luminosity/temperature of galaxy clusters, our prediction can be compared with 
the upcoming cluster redshift catalogues in optical and X-ray bands. Furthermore, we 
discuss the detectability of the biasing and cosmological model dependence from the 
future data. 

Subject headings: cosmology: theory - galaxies: clustering - galaxies: clusters: general - 
galaxies: halos - dark matter - large-scale structure of universe - methods: statistical 



1. INTRODUCTION 

The genus statistics is a quantitative measure of the topological nature of the density field. In 
contrast to the more conventional two-point statistics like the two-point correlation functions and 
the power spectrum, genus is sensitive to the phase information of the density field, and therefore 
plays a complementary role in characterizing the present cosmic structure. 
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Mathematically, genus G(u) is defined as —1/2 times the Euler characteristic of the isodensity 
contour of the density field at the level of v times the rms fluctuations a. In practice this is equal 
to (number of holes) — (number of isolated regions) of the isodensity surface. In cosmology it is 
convenient to use the genus per unit volume, or the genus density g{v), rather than the genus itself. 
The genus density of the random-Gaussian density field is analytically given by 



where o\ = {[V5(x)] 2 ) and a = {[S(x)] 2 ) with () denoting the average over the probability distri- 
bution function (PDF) of the density field (Doroshkevich 1970; Adler 1981; Bardeen et al. 1986; 
Gott, Melott, & Dickinson 1986; Hamilton, Gott, & Weinberg 1986). 

In fact the standard model of structure formation assumes that the primordial density field 
responsible for the current cosmic structure is random-Gaussian. Thus equation (1) is the most 
important analytical result for genus, and has been compared with the observational estimates from 
galaxy and cluster samples in order to probe the random-Gaussianity in the primordial density 
field (e.g., Vogeley et al. 1994; Rhoads et al. 1994; Canavezes et al. 2000). Nevertheless the 
nonlinear gravitational evolution inevitably changes the random-Gaussianity. Matsubara (1994) 
explicitly finds the next-order terms for the genus using the perturbation analysis of the evolution 
of the density field. His formula is in good agreement with the results from N-body simulations 
(Matsubara & Suto 1996; Colley et al. 2000) in a weakly nonlinear regime. The genus in a strongly 
nonlinear regime has been considered using a variety of modelings including the direct numerical 
simulations (Matsubara & Suto 1996), the empirical log-normal PDF for cosmological nonlinear 
density field (Matsubara & Yokoyama 1996), and the Zel'dovich approximation (Seto et al. 1997). 

Those previous studies, however, did not take account of the effect of biasing between luminous 
objects and dark matter. On the other hand, there is growing evidence that the biasing is described 
by a nonlinear and stochastic function of the underlying mass density (Dekel & Lahav 1999). If this 
is the case, the genus for the isodensity surface of luminous objects should be significantly different 
from that of mass. Of course the biasing is supposed to be very dependent on the specific objects, 
and the general prediction is virtually impossible. Therefore in this paper we focus on a nonlinear 
stochastic biasing model of dark matter halos (Taruya & Suto 2000). This model is based on the 
extended Press-Schechter theory (e.g., Bower 1991; Bond et al. 1991) combined with the Mo & 
White (1996) bias formula, and is in reasonable agreement with numerical simulations (Kravtsov Sz 
Klypin 1999; Somerville et al. 2001; Yoshikawa et al. 2001). In this paper, we present an analytic 
expression for genus of dark matter halos by taking account of the nonlinearity in the halo biasing. 
Since recent numerical simulations (Taruya et al. 2001) suggest the stochasticity in halo biasing is 
fairly small on scales ^ 10/i -1 Mpc (see their Fig. 8), we neglect the effect of stochasticity in what 
follows. The one-to-one correspondence between those halos and galaxy clusters is a fairly standard 
assumption, and our predictions for halo genus as a function of halo mass can be translated into 
those of clusters with the corresponding luminosity/temperature. Therefore our predictions are 
useful in future data analysis of the cluster surveys in different bands. 
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The plan of this paper is as follows. In §2, we briefly review the nonlinear stochastic biasing 
model for dark matter halos that we adopt, and derive an analytic expression of genus for halos. In 
§3, we present specific predictions of genus statistics paying particular attention to the dependence 
on the smoothing length, the range of halo mass, cosmological parameters and the redshift. Finally 
§4 is devoted to conclusions and discusses implications of our results. 



2. PREDICTING THE GENUS FOR HALOS 

In this section we derive an analytic expression for genus of halos adopting the nonlinear 
stochastic biasing model (Taruya & Suto 2000). First we briefly review the halo biasing model 
(§2.1). Assuming that PDF of nonlinear mass fluctuations is given by the log-normal distribution, 
we compute the PDF of halo number density fluctuations, -Phaio(^haio) (§2.2). Finally in §2.3 we 
present our analytic expression for genus of halos which properly takes account of both the biasing 
and the nonlinear gravitational evolution for the first time. 



2.1. Halo Biasing Model 

The central idea of the halo biasing model of Taruya & Suto (2000) is to regard the formation 
epoch and mass of halos as the major two hidden (unobservable) variables which lead to the 
nonlinear stochastic behavior of the halo number density field, 5haio> as a function of the underlying 
mass density fluctuation, 5 mass . Applying the extended Press-Schechter theory, Mo & White (1996) 
derived an expression of Jhalo at z smoothed over the scale of R s as a function of 5 mass , the halo 
mass M, and its formation epoch zf. 

£haio = A h ( J R s , z\5 mass ,M, Z{). (2) 

Since the above procedure assumes the spherical collapse model, we adopt the top-hat smoothing 
in computing <5mass and 5haio throughout the paper. Convolving the above expression with the 
Press-Schechter mass function n(M, z; S Cj o) and the halo formation epoch distribution function, 
dp/dzt{zt\M,z) (Lacey & Cole 1993; Kitayama & Suto 1996), Taruya & Suto (2000) derive the 
following conditional PDF of 5 halo for a given S mass , P(S halo \5 mauSS ): 

P(S halo \5 mass ) dd halo = Af- 1 [ [ dM dz f ^-( Zi \M,z) n(M,z;6 Ct0 ), (3) 

J JC(M,Z{) "Zf 

where <5 Ci o is the critical threshold for spherical collapse 3(127r) 3 ' 2 /20 ~ 1.69. The region of the 
integration, C(M,Zf), is given as follows: 

C(M, z { ) = { (M, Zf) | (5 ha i < A h (i? s , z\6 mass , M, z { ) < S haio + d5 hSL i , 

M min < M < M max , z < z { < oo }, (4) 
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where M m \ n and M max denote maximum and minimum of halo mass. The normalization factor N 
is defined as 



Af= dM dz { -^( Zi \M,z)n(M,z;S c ,o). (5) 

JM min Jz OZi 



The joint PDF P(5 mass , <5 ha i ) is simply given by multiplying the PDF of <5 mass , -P(<5 mass ), which we 
assume is log-normal: 

-P(<5mass, <5t, a lo) = ^(^halol^mass) -P(5mass)- (6) 

It is well known that the PDF of the nonlinear mass density field is approximated by the 
log-normal distribution (e.g. Coles & Jones 1991; Kofman et al. 1994). Using this PDF in the 
present context, however, implies that we implicitly assume the one-to-one correspondence of the 
primordial density fluctuation and the evolved nonlinear one (<5 mass ). Recently Kayo et al. (2001) 
show that the scatter around the mean relation between the initial and final density fluctuations is 
significant although the nonlinear PDF is in good agreement with the log-normal. Thus it should 
be noted that the present paper neglects the effect of the scatter or the stochasticity of the mass 
density field. This effect will be examined elsewhere when we compare the predictions of this paper 
with results of N-body simulations (Hikage et al. in preparation). In the same spirit, we below 
do not take into account the stochasticity of Jhalo with respect to <5 mass , and adopt the one-to-one 
mapping between <5haio and <5 mass represented by the mean biasing: 

<5halo = J -P(5halo|<5 m ass) ^halo ^halo = /bias(^mass)- (7) 



In what follows, we use <5haio to indicate 5haio(^mass) = /bias(<5 m ass) in reality. 

In §3.2, we show that our results can be well reproduced from the mean biasing function up 
to the second-order term of <5 mass . In fact, this is important in the sense that our genus expression 
in the next section is applicable to the more general biasing scheme as long as its second-order 
parameterization form is specified. In the present biasing model, we first compute the mean biasing 
numerically according to equation (7), and then fit the data in the range of —0.1 < 5 mass < 0.1 to 
the following quadratic model: 



fbils ^ (^mass) — ^l^mass + ^ (^mass "mm ) (8) 



which satisfy the relation (f^^\5 mass )} = (a mm is defined in eq.[10] below). The fitted values 
for the linear and the second-order biasing coefficients, 61 and 62, are listed in Table 1 for several 
models. 

Figure 1 plots /bias(^mass) and f^g\5 mass ) in solid and dotted lines, for different M m - m , R s , and 
z. Since we are primarily interested in the galaxy cluster scales, we consider the following values 
of the parameters; the minimum mass of halo M m j n = 1O 13 /i _1 M and W lA h~ 1 M Q , the smoothing 
length i? s = 30/i _1 Mpc and 50/i~ 1 Mpc, z = and 1 (h denotes the Hubble constant in units of 
100 km-s _1 -Mpc _1 ). While we set the maximum halo mass M max = W 16 h" 1 M & , the larger value 
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of M max does not change the results below since the number density of such massive halos is quite 
small. Our predictions below are computed in two representative Cold Dark Matter (CDM) models; 
Lambda CDM (LCDM) with (Q , A , a 8 , h) = (0.3,0.7,1.0,0.7), and Standard CDM (SCDM) with 
(J7o) -V)) 08> ft,)=(1.0,0.0,0.6,0.5). The density parameter Qq, dimensionless cosmological constant, 
Ao, and the top-hat mass fluctuation at 8/i _1 Mpc, erg are normalized according to the cluster 
abundance (Kitayama & Suto 1997). 

Figure 1 clearly indicates the degree of nonlinearity in our biasing model. For the range of 
parameter values of our interest, the quadratic fit provides a reasonable approximation to 

the mean biasing. As M m \ D and z increase, nonlinearity in the biasing becomes stronger (see 62 in 
Table 1 ), and one expects the more significant departure of the resulting genus density from the 
random-Gaussian prediction (eq.[l]) as we will show below. 



dS mass ■ 



2.2. Probability Distribution Function of Nonlinear Mass and Halo Density Fields 

PDF of the objects of interest is the most important and basic statistics in computing their 
genus. The PDF of halos in our biasing model is derived in this subsection. 

Standard models of structure formation inspired by the inflation picture assume that the 
PDF of the primordial density field obeys the random-Gaussian statistics. While the nonlinear 
gravitational evolution of the density field distorts the primordial random- Gaussianity, the log- 
normal PDF is known to be an empirical good approximation to the nonlinear mass density field 
<5 m ass (see discussion in the previous subsection): 

p (A W e _ 1 f {m[(l + W)yi + <Tmm 2 ]} 2 " 

^rnassVOmassJOOmass — r- — ; — j- ^ 6Xp „ 

"mass 

) v/Z7T m( 1 + @ mm ) L V ^ ' ^mm ) 

(9) 

The variance a mm (R s , z) of the mass fluctuations at z smoothed over the scale R s is computed 
from the nonlinear power spectrum P(k, z) convolved with the specific window function W 2 (kR s ): 

/dh 
— k 2 P m (k,z)W 2 (kR s ) . (10) 

In what follows, we adopt the fitting formula of the nonlinear CDM power spectrum P m (£;, z) 
(Peacock & Dodds 1996). Throughout the paper, we adopt the top-hat window function: 

3 

W(x) = ^-(sinx — xcosx). (11) 
x 6 

The resulting PDF in this prescription proves to be in excellent agreement with N-body simulation 
results (Kayo et al. 2001). 

One interpretation of equation (9) is that the nonlinear mass density field 5 ma ss has the following 
one-to-one mapping from a random-Gaussian field 5 RG with unit variance: 

<5mass = /iii[<5rg] = 1 ,> exp [^/M 1 + CW^rc] ~ 1- (12) 



_ 2 

(J mm 



Vi + ct ; 



1 

mm 
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Physically speaking, the above 5 RG corresponds to the primordial density field apart from the 
arbitrary normalization factor. Of course this mapping cannot be strict since (5 mass should not be 
determined locally. Nevertheless this interpretation provides an interesting and useful method to 
compute genus of nonlinear mass density field as discussed by Matsubara & Yokoyama (1996). 

We extend this idea in our biasing model. Since our halo density field is given by the mean 
biasing (7) and we adopt the log-normal PDF (9) for <5 mass , we can easily obtain the PDF for 

<5halo = /bias(<5mass) = /bias (fin (<W) as follows: 



-Phalo(<5halo)<#>halo = -Pmass[/ b 7as( ( ^halo)] 



mass U bias V"iia,io/j j r I r 
u J bias yJm 

1 



^halo (13) 

<5mass —/bias (^halo ) 



x exp 



/bias(/b7as(^alo))(l + /^(W)AMl + 

{ln[(l + / b li(5haio))v / T+^7]} 



21n(l + <7 mm 2 ) 



^halo- (14) 



In the above, / bias (5 m ass) denotes the derivative respect to <5 m ass, and / bi * s is the inverse of the 
monotonic function /bias- 

Figure 2 plots the PDF for dark halos in solid lines, respectively, employing the same sets of 
parameters as Figure 1. In order to separate the changes of the PDF due to the nonlinear biasing 
and the nonlinear gravitational evolution, we compare the dark halo PDF with the log- normal PDF 
(eq:[9]; dotted lines). Also for that purpose the variance cr mm of the log-normal in equation (9) is 
adjusted to the corresponding cihh of the halos: 

o"hh 2 = J J <5haio 2 -P(<5 mass i C^halo, (15) 

in plotting the log-normal PDF curves, Note that the quantity a na is somewhat different from 
variance of the averaged biasing, ([5haio(<5mass)] 2 ) and contains the nonlinearity and the stochasticity 
of biasing. 2 

Figure 2 shows the sensitivity of PDFs for the halo mass M m i n and the cosmological parameters, 
especially for the fluctuation amplitude as- The dark halo PDF significantly differs from the log- 
normal PDF even using the same value of variance, which clearly indicates the importance of the 
nonlinearity in the halo biasing. In some models, the PDF at 5 nai \ = — 1 diverges since <5h a io i n our 
model reaches —1 even though <5 maS s is greater than —1. Then /bi as is not a one-to-one mapping 
between 5 na i and <5 mass and the expression (13) becomes singular there. Of course, this divergence 
occurs only at 5h a io = — 1 5 an d does not affect our results for 5haio > _ 1- 



"According to Taruya & Suto (2000), quantities Uhh and ([5haio(<5mass)] ) are related in the following way 



{[<5halo(<5mass)] ) = Ohh 



1 + 

i + e + 4 



nl ' c scatt 

where e n i and e sca tt respectively denote the degree of nonlinearity and stochasticity. 
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2.3. Genus Curve For Dark Halos 



As we describe in the previous subsection, our halo biasing model results in the one-to-one 
mapping between S nG and <5 ha io = /bias(^mass) = /bias (/invito))- Then using the prescription of 
Matsubara & Yokoyama (1996), we can derive the genus density explicitly. 

For the nonlinear mass density <5 mass with the log- normal PDF (9), the genus density is com- 
puted as (Matsubara & Yokoyama 1996): 



5mass( 1/ m) — 9m AX 6Xp 



{ln[(l + V m (J mm )\/l + Q "mm 2 ]} 

21n(l + cr mm 2 ) 

21 \2 



21 



1 _ {ln[(l + i/ m q- mm )\/l + o-mm 2 ]} 
ln(l + <J mm 2 ) 



(16) 



where the level of isodensity contour u m is given by ^mass/^mm- The quantity g M Ax is the maximum 
value of the genus: 



9max — 5 f mass( < ^MAx/°"mm) — 



l,m 



(2vr) 2 [3(l+a mm 2 )ln(l + a mm 2 )] 3 / 2 ' 
at <5 MAX = 1/ \/l + o"mm 2 - 1 with 

/rlh 
^k 4 P nl (k,z)W 2 (kR s ) . 



(17) 



(18) 



We transform this genus density using the mean biasing relation <5haio = /bias(<5 m ass)- Then, at 
a given threshold = 5haio/chh> our final expression for the genus density of dark halos becomes 
as follows: 



S'haio(z'h) = 9max exp 



21n(l + a mm 2 ) 
ln(l + cr mm 2 ) 



(19) 



with cr hh defined by (15). Here, g MAX = ^ a io[/bias(<Wx)/o-hh] = 5mass(<Wx/o"mm) is given also by 
equation (17) and is listed in Table 2 for some models. 



3. PREDICTIONS OF GENUS FOR HALOS IN CDM MODELS 

In this section we present several model predictions of the genus density for dark matter halos 
with particular emphasis on its parameter dependence. Then we examine to which extent our 
predictions are reproduced with a perturbation model on the nonlinear gravitational growth and 
the biasing. 
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3.1. Parameter Dependence of the Halo Genus 

Consider the halo mass dependence first. As we argued in §2.1, the maximum mass does not 
affect the results as long as M max > 1O 16 M . They are, however, very sensitive to the minimum 
mass M m j n . This mass dependence is shown in Figures 3 and 4 in LCDM and SCDM models, 
respectively, where we plot the normalized genus density (i.e., g(y) in units of its maximum value 
5max)- As M m i n increases, the nonlinearity of both the halo biasing and the mass fluctuation 
amplitude becomes stronger (see Fig.l and Table 1) which results in the significant departure from 
the random-Gaussian prediction (dashed lines). 

Also the genus is fairly dependent on the smoothing length R s (Figs 3 and 4, and Table 2). 
While the genus with R s ^ 100/i -1 Mpc approaches the random-Gaussian prediction, the deviation 
is detectable with smaller R s . We discuss the detectability of the signature in the last section. 

In contrast, the redshift dependence is rather weak; neither the amplitude and the shape of 
the genus density evolves much, at least between z = and 1 (which is the relevant range of the 
redshift for galaxy clusters). This is interpreted as a result of the two competing effects; the stronger 
clustering due to the halo biasing at higher z and the weaker amplitude of the mass fluctuations (see 
also Taruya Sz Yamamoto 2001). 

The above features apply both to LCDM and SCDM. The comparison of Figures 3 and 4 
indicates that the predictions in SCDM model is slightly closer to the random-Gaussian. This is 
simply because these models are normalized according to the cluster abundance and thus the mass 
fluctuation amplitude at z = is smaller in SCDM (erg = 0.6) than in LCDM (og = 1.0). 



3.2. Comparison with the Perturbation Analysis 

So far we have presented our predictions using the fully nonlinear model of both the halo 
biasing and the mass fluctuations. As long as galaxy clusters are concerned, however, we are 
mainly interested in fairly large R s . Therefore the perturbation analysis may be a reasonable 
approximation. This line of studies is pioneered by Matsubara (1994) in the weakly nonlinear 
evolution of mass density field, and also by Matsubara (2000) including the effect of biasing. We 
will perform the perturbation analysis of our full model predictions. This attempt is instructive 
in understanding the origin of the departure from the random-Gaussianity. Also the resulting 
perturbation formula becomes useful since it is readily applicable to a wide range of the biasing 
models, in addition to our halo biasing, once their second-order biasing coefficients are given. 

Matsubara (1994) derives the analytical formula for the genus of the mass density field up to 
the lowest-order term in the (linear) mass variance a : 



9{») 



1 



(27T) 



Wo-) 



y ) { Hi{v) + 



S ( ) 



(20) 



-9- 



where represents the skewness parameters denned by 

W 3 ) 



S<°> = (21) 
a 



= 3 



4 ff? 



(22) 



5(2) = _9((w 1 wxvM)^_ (23) 

We expand equation (16) with respect to cr mm and find that it reduces to 

5mass(^m) = S^RG^m) | 1 + ^H 3 (u m )a mm + C((7 mm 2 ) 1 , (24) 



where 



1 \ 3 /2 2 



SRO („ m) = _^ (l-i)- («) 

This implies that all the above skewness parameters are equal to 3 in the log-normal PDF. 

If we further include the halo biasing effect using the second-order fit <5h a io = f^^\S mass ) given 
by equation (8) , we obtain the perturbation formula of our halo genus up to <7hh as 

ffhaloK) = Shalo.oK) jl + ^ ~^~ ff 3 ^ V ^ h ^j a hh + C((7 hh 2 )| , (26) 



2 \ f / \ 2 i 

"var^h \ II-, / "var^h 



where 

^ j iimif-mi' (27> 

with 6 var = cr hh /c7 mm (Table 1). 

If the biasing is linear and deterministic, 6 V ar = b\ and the skewness parameters reduce to 



sP-s£"-flP-£(l + £], (28) 



which agrees with equation [4.82] of Matsubara (2000). 

Figure 5 compares those perturbation results with our full model predictions; Top panel shows 
the normalized genus in the log- normal model (solid lines; eq.[16]) and in the perturbation model 
(dotted lines; eq.[24]). Dotted lines in Middle panel refer to the results combining the weakly non- 
linear evolution model (eq.[24]) and the second-order biasing model (eq.[8]), while in Bottom panel 
they refer to the results combining the log-normal mass PDF (eq.[16]) and the second-order biasing 
model (eq.[8]). Those are to be compared with our full predictions (solid lines). In the perturbation 
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results (eqs.[24], [26]) we define g M Ax respectively as the maximum value in the set of the plotted 
genus respectively. This comparison indicates that the second-order biasing model provides a very 
good approximation for the parameters of interest, also that the lowest-order correction to the 
nonlinear mass evolution dominates for large smoothing lengths R s ^ 50/i -1 Mpc. Therefore we 
conclude that our perturbation formula (26) is in practice a useful and reliable approximation to 
the genus of galaxy halos, which is applicable even beyond our particular biasing model discussed 
here. 



4. SUMMARY AND DISCUSSION 

In the present paper, we have presented an analytic model of genus for dark matter halos 
adopting the stochastic nonlinear halo biasing model by Taruya & Suto (2000), assuming that the 
primordial mass density field obeys the random-Gaussian statistics. This is the first attempt to 
predict the genus statistics simultaneously taking account of the nonlinear nature in the biasing 
and of the nonlinear gravitational evolution of the underlying mass distribution. In the remainder 
of this section, we discuss several implications of our model predictions. 

First of all, we note that our model applies only for dark matter halos identified according to 
the spherical collapse model strictly speaking. Nevertheless this picture is now widely accepted as 
an empirical model for galaxy clusters in the universe. In fact, the excellent agreement between 
the predicted halo abundance and the observed cluster abundance (e.g., Kitayama & Suto 1997; 
Suto 2000 and references therein) justifies the empirical one-to-one correspondence between the 
theoretical halo and the observed cluster in a statistical sense. Adopting this conventional view, 
our predictions for dark matter halos can be readily translated into those for clusters on the basis 
of the mass-temperature and mass-luminosity relations at the present epoch, for instance, as 

k - T = ls (it*) 1/3 do4wr nfkev> <29> 

and 

^ - (iotoF"" 1042 ^ sec_i ' <3o> 

where k B is Boltzmann constant and A c is the mean density of a virialized cluster (Kitayama & 
Suto 1996; Suto et al. 2000). Therefore our model predictions can be tested against future cluster 
observations in various wavebands. 

Next it should be stressed that our current model is based on the one-to-one mapping of the 
primordial density field that obeys the random-Gaussian statistics. In this approximation, if we 
express our genus predictions as a function of V defined through the volume fraction /: 
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they simply reduce to the random-Gaussian prediction (eq.[l]). Most previous studies used this V in 
order to remove the (unknown) distortion of the PDF due to the gravitational nonlinear evolution. 
Nevertheless we prefer to adopt the correct definition of v = 5/a because the departure from 
the Gaussian prediction (eq. [1] ) provides also an important measure of the gravitational nonlinear 
evolution and the nonlinear biasing in more advanced approach than what we took here. For 
instance, the perturbative results in weakly nonlinear regime (Matsubara 1994) do not reduce to 
equation (1) even in terms of 17, and also the effect of stochasticity in biasing, which is ignored 
in the present analysis, would add another feature which violates the simple scaling idea behind 
the use of v. As we have shown in the present paper, the genus as a function v can be a useful 
probe of the nonlinear gravitational evolution and the biasing of objects and mass, as long as the 
random-Gaussianity of the primordial density field is correct. The comparison between g(p) and 
g{v) is discussed in detail by Colley et al. (2000). 

In our genus model, the overall amplitude of genus for dark halos g MA x depends only on 
mass density spectrum. Therefore the observed amplitude of the cluster genus may provide direct 
information on the mass fluctuation in principle. However this fact is heavily dependent on our 
assumption of one-to-one mapping of density fluctuation in equation (7), and needs to be checked 
with future numerical simulations. 

Finally we briefly discuss the future detectability of the departure of the halo genus from the 
random-Gaussian prediction. We use the top-hat smoothing R s = 50/i -1 Mpc and compute our 
model predictions for M min = 1O 13 /i -1 M and 1O 14 /i _1 M in SCDM and LCDM models. Figure 
6 plots the total number of genus G(u); Upper and Lower panels are normalized so as to roughly 
correspond to the Sloan Digital Sky Survey (SDSS) spectroscopic and photometric cluster samples, 
respectively. We do not take into account the light cone effect (Yamamoto & Suto 1999; Suto et al. 
1999) but rather show the results at the mean redshift z = 0.1 and z = 0.3. Taruya & Yamamoto 
(2001) find that the light-cone effect is not so significant for genus statistics up to z ~ 1. The 
SDSS is expected to detect N spect ~ 1000 clusters from spectroscopic data up to z ~ 0.2 and 
-Wphoto ~ 5000 from photometric data up to z ~ 1. Assuming that errors of genus are dominated 
by the sampling noise, we assign the error independently of v as 5G = \/g M AxVcS, where V e fi is the 
effective volume of the sample. We further assume that V e s is proportional to the number of clusters. 
If we adopt the value for the Abell cluster sample of Rhoads et al.(1994), ~ 250 Abell clusters in 
V c ff = 9.5 x 10 6 (/i _1 Mpc) 3 , V e s for SDSS spectroscopic and photometric samples is estimated to 
be 4 and 20 times larger than that of Rhoads et al.(1994), respectively. The amplitude of the error 
bars estimated in this way indicates that the departure from the random-Gaussianity should be 
indeed detectable especially in photometric redshift data and that one can even discriminate among 
some model predictions with the upcoming data of galaxy clusters. 

We plan to improve our model predictions by including the stochasticity in the biasing model, 
the light-cone effect, and redshift-space distortion. Also we are currently working on evaluating the 
genus of halos identified from the cosmological N-body simulations. These results and comparison 
will be reported elsewhere. 
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Table 1. The linear and second-order biasing coefficients, b\ and &2(eq.[8]), respectively, and the 
ratio of the rms fluctuations between dark halos and mass 6 V ar (= ^hh/^mm) for different 

cosmological models, z, M m i n , and i? s . 



Model 


z 


M min [h- l M Q ] 


tfs^Mpc] 


h 


b 2 


^var 


LCDM 





10 13 


30 


1.77 


0.23 


2.12 


LCDM 





10 13 


50 


1.76 


0.20 


1.87 


LCDM 





10 13 


100 


1.76 


0.14 


1.84 


LCDM 





10 14 


30 


2.72 


3.26 


4.90 


LCDM 





10 14 


50 


2.92 


3.20 


3.34 


LCDM 





10 14 


100 


3.00 


3.70 


3.29 


SCDM 





10 13 


30 


1.94 


0.51 


2.07 


SCDM 





10 13 


50 


1.94 


0.77 


2.04 


LCDM 


1 


10 13 


30 


2.85 


2.78 


3.22 


LCDM 


1 


10 13 


50 


2.90 


2.60 


3.08 



Table 2. The maximum values of the genus g M Ax for different cosmological models, z, and R s . 



Model 


z 


Rslh^Mpc] 


5max [(100/i- 1 M P c)- 3 ] 


LCDM 





30 


4.75 


LCDM 





50 


1.49 


LCDM 





100 


0.38 


LCDM 


1 


30 


3.56 


LCDM 


1 


50 


1.13 


LCDM 


1 


100 


0.29 


SCDM 





30 


9.47 


SCDM 





50 


3.83 


SCDM 





100 


1.48 


SCDM 


1 


30 


7.66 


SCDM 


1 


50 


3.12 


SCDM 


1 


100 


1.22 
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30 
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50 
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LCDM 
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LCDM 

z = 1 

M„ ; =10 13 h-W 




30 

Fit to eq.[8] _ 

50 

Fit to eq.[8] 
I M ill 



flglh^Mpc] 
30 

Fit to eq.[8] _ 

- 50 

- Fit to eq.[8] 



j i 



0.1 



1+5. 



0.1 



l + <5. 



10 



Fig. 1. — The mean biasing function /bias(^mass)(soZic?) and its fitting up to the second-order of 
<5mass> /bias^^massX dotted). The red and blue lines indicate the results with R s = 30/i _1 Mpc and 
50fc _1 Mpc. Upper-left: LCDM model, z = 0, and Af min = 10 13 /» -1 M o ; Upper-right: LCDM 
model, 2 = 0, and M min = lO 14 /^ 1 ^/© ; Lower-left: SCDM model, z = 0, and M min = lO^/i^M© 
; Lower-right: LCDM model, z = 1, and M min = lO^/i^M©. We set M max = lO 16 /i _1 M . 
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LCDM M . =10 13 h- 1 M o 

min fc> 

z = R s [h- l Upc][ a hh ] 
- 30 [ 0.67 ] 

log-normal 

50 [ 0.33 ] 
log-normal 




SCDM 

z =0". 




M min =W^h-W Q 
R s [h->Upc][ a hh ] 

30 [ 0.28 ] 

log-normal 
50 [0.14] 

log-normal 



j i i ± 



— i — i — i — | — i — i — i — | — i — i — i — | — i — i — r 

LCDM M.=10 u h~ l M n 

mm Cy 

z = /? s [/i- 1 Mpc][ a hh ] 

30 [ 1 .54 ] 

log-normal 
50 [ 0.60 ] 
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M min =lOi 3 h-iM 
* s [h-iMpc][ a hh ] 

30 [ 0.63 ] 

log-normal 

— 50 [ 0.34 ] 
log-normal 







1 

6 







1 
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Fig. 2. — Probability distribution functions for dark halos with R s = 30h Mpc and 50h Mpc. For 
reference, dotted lines indicate the log-normal distribution with the variance of the corresponding 
a hh . Upper-left: LCDM model, z = 0, and M min = 10 13 h,- l M & ; Upper-right: LCDM model, z = 0, 
and M min = 10 14 h- 1 M Q ; Lower-left: SCDM model, z = 0, and M min = lO 13 /i _1 M ; Lower-right: 
LCDM model, z = 1, and M min = W^h^MQ. In all panels M max = lO 16 /i _1 M . 
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Fig. 3. — Normalized genus density 5 f (^)/fl , MAx 

for dark halos in LCDM model with R s = 30/i _1 Mpc 
(red), 50/i _1 Mpc (blue), and 100/i _1 Mpc (green). For comparison, Upper-left panel plots the nor- 
malized genus density for dark matter. The other panels assume M max = 10 h M@. Upper-right: 
z = and M min = 1O 15 /i _1 M ; Middle-left: z = and M min = 10 14 /i _1 M o ; Middle-right: z = 
and M min = lO^/i^M©; Lower-left: z = 1 and M min = W u h- 1 M Q ; Lower-right: z = 1 and 
M min = lO 13 /*- 1 ^^. 
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Fig. 4. — Same as Figure 3 for SCDM model. 
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Fig. 5. — Comparison of <?(^)/(7max between the full and perturbation models. The predictions 
assume z = in LCDM model, M min = 10 14 /i~ 1 M Q , M max = 10 16 h~ 1 M Q , and three different 
smoothing scales; R s = 30/i _1 Mpc (red), 50/i _1 Mpc (blue), and 100/i _1 Mpc (green). Top: mass 
in the log- normal model (solid; eq.[16]) and in the perturbation model (dotted; eq.[24]). Middle: 
dark halos (solid; eq.[19]) and the perturbation formula of our halo genus with the second-order 
biasing model (dotted; eq.[26]). Bottom: dark halos (solid; eq.[19]) and the results combining the 
log-normal mass PDF and the second-order biasing model(dotted; eqs.[8], [19]). 
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Fig. 6. — Expected total number of genus for surveys corresponding to N ~ 1000 ( Upper panel) and 
N ~ 5000 (Lower panel); M min = lO 13 /i -1 M (solid) and M min = lO 14 fr -1 M (dotted) in LCDM 
(red) and SCDM (blue) models. All the curves adopt top-hat smoothing of R s = 50/i~ 1 Mpc. 



